Statistical-mechanical theory of ultrasonic absorption in molecular liquids 



o 
o 

O 

Q 



o 

£ 

o3 



G 

O 

o 



> 
o 

(N 

O 

-i— > 



i 

C 

o 
o 



X 



Department of Theoretical Study, 



Alexander E. Kobryn 
Institute for Molecular Science, Myodaiji, 



Okazaki, Aichi 4-44-8585, Japan 



Fumio Hirata 

Department of Theoretical Study, Institute for Molecular Science, 
and Department of Functional Molecular Science, 
The Graduate University for Advanced Studies (SOKENDAI), Myodaiji, Okazaki, Aichi 444-8585, Japan 

(Dated: October 10, 2006) 

We present results of theoretical description of ultrasonic phenomena in molecular liquids. In par- 
ticular, we are interested in the development of microscopical, i.e., statistical-mechanical framework 
capable to explain the long living puzzle of the excess ultrasonic absorption in liquids. Typically, ul- 
trasonic wave in a liquid can be generated by applying the periodically alternating external pressure 
with the angular frequency that corresponds to the ultrasound. If the perturbation introduced by 
such process is weak - its statistical-mechanical treatment can be done with the use of the linear re- 
sponse theory. We treat the liquid as a system of interacting sites, so that all the response/aftereffect 
functions as well as the energy dissipation and generalized (wave-vector and frequency dependent) 
ultrasonic absorption coefficient are obtained in terms of familiar site-site static and time correla- 
tion functions such as static structure factors or intermediate scattering functions. To express the 
site-site intermediate scattering functions we refer to the site-site memory equations in the mode- 
coupling approximation for the first-order memory kernels, while equilibrium properties such as 
site-site static structure factors, direct and total correlation functions are deduced from the integral 
equation theory of molecular liquids known as RISM or one of its generalizations. All the formalism 
is phrased in a general manner, hence the obtained results are expected to work for arbitrary type 
of molecular liquid including simple, ionic, polar, and non-polar liquids. 

PACS numbers: 05.20.Jj - Statistical Mechanics of Classical Fluids, 61. 25. Em — Molecular Liquids 
Keywords: Ultrasonic Absorption Coefficient, Integral Equation Theory, Memory Equation, 



I. INTRODUCTION 

Sound velocity and absorption are important probes 
of physical and chemical processes taking place in liquid 
and liquid mixtures [1, 2]. Sound velocity reflects density 
fluctuations occurring in solutions, thereby, it is closely 
related to the compressibility, which is a response func- 
tion of the density to the mechanical perturbation, or 
pressure. The property is probably the best measure for 
the mechanical stability of solutions. On the other hand, 
sound absorption probes the energy dissipation caused 
by a variety of irreversible processes taking place in solu- 
tions: for example, structural relaxation of liquids, con- 
formational transition of molecules in solution, and so 
forth. No wonder why people have employed the method 
for long time to investigate the structure and dynamics 
in solutions. There has been enormous amount of ex- 
perimental data accumulated in the filed of science and 
technology. More importantly, the method seems to be 
finding new horizon in different fields of scientific research 
[3-7], not mentioning about the usage in medicine and 
pharmacy [8-17]. Especially important among the appli- 
cations is that to the stability and conformational transi- 
tion of biomolecules such as protein. The sound velocity 
and the adiabatic compressibility have already been on 
duty to clarifying the mechanical stability of protein in 
solution [18-21]. At the frontier in the science, sound 
velocity has been employed to identify the conforma- 
tional transition from native to denatured states of pro- 



tein (for the effect of ultrasound on biological molecules 
and macromolecules see, e.g., Refs. [22-35]). It won't be 
so long before the sound absorption associated with the 
conformational relaxation is determined by the acoustic 
measurement. The method seems appropriate to investi- 
gate the large scale fluctuation or the collective dynam- 
ics characteristic to the macromolecule, because the fre- 
quency range of such motion is covered by that of sound 
wave. However, there exists a high barrier for the method 
to be overcome in order to be applied to the field. That 
is the method or theory to analyze the experimental ob- 
servations. 

Traditional ways of interpreting data from acoustic 
measurements are based on the hydrodynamic and/or the 
phenomenological theories of relaxation [36-45] exempli- 
fied by that of Debye-type and of stretched-exponential- 
type. Those theories inevitably require phenomenological 
and/or molecular models for the relaxation, which some- 
times are far from reality. One of the most successful 
theories of the sound absorption is that for water, pro- 
posed by Hall more than half a century ago [46, 47]. Hall 
has assumed that water is an equilibrium mixture of two 
states: monomers having densely packed structure and 
hydrogen-bonded clusters which have an ice-like open- 
packed structure. The equilibrium will shift toward the 
monomers by acoustic pressure due to the difference in 
the molar volume between the two components of the 
mixture. The energy dissipation associated with the re- 
laxation process is probed by the sound absorption. Hall 
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applied the transition state theory to the two-state model 
in order to explain the sound absorption by liquid water. 
The treatment should be remarked as the most important 
contemporary achievement in theories of sound absorp- 
tion by liquids, because it is the first theory to relate the 
acoustic process with the microscopic model of the liquid 
structure. 

Unfortunately, the progress seems stopped at the level 
of Hall as far as the microscopic treatment of the sound 
absorption is concerned with. The reason is simply be- 
cause Hall's theory employs a structural model for the liq- 
uid. It means that the theory requires a model for liquid 
structure for each liquid or solution, and phenomenolog- 
ical or empirical parameters associated with the model, 
such as the molar volume. The difficulty of construct- 
ing a model for liquid structure will be easily understood 
by taking the water-alcohol mixture as an example. All 
those are suggesting that the theory of sound absorption 
is necessary to start from a Hamiltonian level or a molec- 
ular interaction, not from structural models of liquid. 

One may think it will then be appropriate to employ 
the molecular dynamics simulation, because the method 
starts definitely from the Hamiltonian model. However, 
the problem is formidably difficult for the method, be- 
cause it touches both the hydrodynamic and thermody- 
namic limits. Just imagine the typical wave length of 
the ultrasonic wave, and how many molecules should be 
involved in a wave packet. But if the simulation is tech- 
nically feasible, it does not help much, because the anal- 
ysis of trajectory will require the statistical mechanics of 
molecular liquids. Without having a method of analysis, 
the trajectory is just a garbage. 

The present paper aims to construct a molecular the- 
ory of sound absorption based on the latest development 
in the statistical mechanics, or the site-site representa- 
tion of the generalized Langevin theory. The theory 
combines two elements in the theoretical physics: the 
reference interaction site model (RISM) and the general- 
ized Langevin equation (GLE). The RISM theory is an 
integral equation method to describe equilibrium struc- 
ture of molecular liquids in terms of the site-site pair 
correlation functions, from which all the thermodynamic 
quantities can be derived [48-51]. The theory and its 
three dimensional generalization have been applied suc- 
cessfully to almost entire spectrum of chemical and phys- 
ical processes taking place in solution from chemical re- 
actions to the molecular recognition by protein [51-53]. 
Among the applications, the results for the isothermal 
compressibility of the water-butanol mixture [54] are of 
particular interest in conjunction with the present top- 
ics, since it is closely related to the mechanical stability 
of solutions. The theory could have successfully repro- 
duced the concentration and temperature dependence of 
the compressibility of the solution observed experimen- 
tally, and it could have been able to describe the behavior 
in terms of the liquid structure. In that sense, the theory 
has already demonstrated its capability to explain the 
molecular mechanism of sound velocity. However, it is 



not complete for describing the entire physics of sound 
propagation, because it is in general accompanied by en- 
ergy dissipation, or attenuation of sound wave. This is 
the reason why we employ the GLE theory - a statistical 
mechanics theory for irreversible processes. In past two 
decades, we have developed a method which combines 
RISM and GLE [51]. The method could have described 
a variety of irreversible processes with great success, in- 
cluding the ion transport [55], dynamics of water [56, 57], 
viscosity [58], dielectric relaxation [59], translational and 
rotational dynamics of a molecule in solutions [60-62]. 
Especially important in relation with the present topic 
is the contribution made by Yamaguchi et al. [63] who 
formulated a theory for ultrasonic vibration potential, or 
the coupling of acoustic and electrostatic perturbations 
in polar liquids. This work provides a good guide for 
the development of a theory of ultrasonic phenomena in 
liquids. 

In the present paper, we propose a theory for sound 
absorption in molecular liquids based on combination of 
the RISM theory, the GLE theory, and the linear re- 
sponse theory. Starting from definition of the pertur- 
bation Hamiltonian which describes the coupling of the 
sound wave and a liquid system, we derive an expression 
for the sound absorption in the most general form. In 
the choice of parameters for abbreviated description we 
recognize that the traditionally employed set of fluctua- 
tions of local site densities and longitudinal site currents 
is not enough to serve our purposes in this case and has 
to be extended by adding fluctuations of the total energy 
density. The latter is really necessary if one wants to 
obtain in the expression for ultrasonic absorption coef- 
ficient terms related to the thermal conductivity of the 
system. We then apply the theory to the simplest case 
of a liquid of spherical molecules and derive expressions 
in the hydrodynamic limit in order to make contact with 
the result from classical hydrodynamics. 

The paper is organized as follows. In Section II we spe- 
cify the system, define basic dynamical quantities that 
we use for its description, and write their equations of 
motion. In Section III we provide derivation of the per- 
turbation Hamiltonian for the case of not too strong ul- 
trasonic vibrations so that the result of derivation can 
be used in the theory of linear response. Application of 
the theory of linear response to our system is described 
in Section IV, where we calculate linear responses of lo- 
cal particle number, charge and mass densities for the 
perturbation caused by the propagating ultrasonic wave. 
All responses are written in terms of familiar static and 
time correlation functions. A major step to calculate 
ultrasonic absorption coefficient is made in Section V, 
where we first derive an expression for the energy dis- 
sipation, and then use it to define the generalized, i.e., 
wave-vector and frequency dependent ultrasonic absorp- 
tion coefficient in a liquid. For that purpose one may 
require the knowledge of site-site intermediate scattering 
functions. Therefore another and rather big portion of 
the material is dedicated to the problem of their calcu- 
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lation and is put in Sections VI and VII. The remaining 
part of the paper is about detailed analysis of our find- 
ings. In particular, we consider the hydrodynamic limit 
of our result for the case of one-component simple liquid 
in order to compare it with the one from the continuous 
media models [1, 2, 36]. The required steps and their 
consequences are described in Section VIII. Finally, ad- 
vances, deficiencies, remedies and open questions of our 
treatment are summarized in Section IX. Auxiliary ma- 
terial is put into Appendices. 



II. DEFINITION OF THE SYSTEM AND 
VARIABLES 

Let us give first generic definitions of major quantities 
for a mixture of A'f-component classical molecular liquid 
confined in volume V. We suppose that each component 
consists of particles (and indices p, v = 1, . . . , A4 are 
used to label molecular species), hence the total num- 
ber of particles is — N. According to this, x M = 
N^/N should be particle relative concentrations (or mo- 
lar fractions). Each molecule consists of A M sites (could 
be atoms) having charges and masses mJJ (and indices 
a, 7 = 1, . . . , Ap are used to label molecular sites). Time- 
dependent positions and momenta of individual sites will 
be denoted then as r™i(i) and p^ A (t) = Tn^v^^i), re- 
spectively, with v"j(i) being site velocities. 

Let us suppose that the unperturbed system in con- 
cern is described by the Hamiltonian H, and that the 
perturbation introduced by the external pressure P is 
H'. Microscopic expression for H can be written as the 
sum of the kinetic ffkm an d interaction H lnt parts as 

H = H kin + H int , (2.1a) 

M N„ A„ r /.y|2 

^=1 i= l a=l 
M A M M N v A„ 

H ** = EEEEEE C (KliM 

fi—l i—l a — l v—1 j—1 7—1 



dH 



r M,iW- r 3jW. and ^,ij(l r ^,ij(*)l) is 



where r^-^) 

potential of interaction between sites, also it is assumed 
implicitly that if v = p,, then j ^ i (in the following, this 
condition will be indicated by the prime mark next to the 
summation symbol). At present, a detailed knowledge of 
structure for 4>'^'l,ij(\ r< ^2,ij(^)\) ^ s not required. For sim- 
plicity of notations we will denote the triplicate sum by 
a symbol as follows 



M n„ A M 

EEE 

\l— 1 i—l a—l 



E< 

® 



M N v A„ 

EEE 

v—l j — 1 7—1 



E 



(2.2) 



Equations of motion for individual site positions and mo- 
menta are the Hamilton equations 



dH 



KM 



(2.3b) 



They are used to derive microscopic equations of motion 
for dynamical variables. Microscopic single-particle den- 
sity and momentum density are defined, respectively, as 

Pl t {v-t) = 5 (r-r« ,(*:)), (2.4a) 
P^M) = P°i(*)5(r-r° (2.4b) 

with 5(r) being the 3d Dirac <5-function. Next we intro- 
duce collective variables such as particle number-density 
p(r;t), particle charge-density p c (r;t) and particle mass- 
density p M (r;t), for which we will use common designa- 
tion p(r; t), so that 

P(r;i) = {p(r;i),p c (r;i),p M (r;i)} 

= ^^(r;i){l,^,m«} 5 

® 

= 5>S,i(r;f)T°, (2.5) 

si 

where is used to denote the set {1, z",m"}. For the 
system in equilibrium (p(r;t)) — p, (p c (r;t)) — and 
(p M (r;t)) = p M , where p is the mean value of the parti- 
cle number-density and p M is unperturbed mass-density. 
Angular brackets (...) denote an appropriate statistical 
average: for example, the Gibbs canonical ensemble ave- 
rage. 

Similarly, the particle number-density current j(r;i), 
charge-density current j c (r;t), and mass-density current 
j M (r;t), for which we also will use one common designa- 
tion jj(r; t), are defined as 

J(r;t) = {j(r;t),jo(r;t),j M (r;t)} 

= Y,Pt,^)^. (2.6) 

® m 

It is worth to note, that j M (r;t) is nothing but the 
momentum-density, i.e., j M (r;£) = p(r;t). Its spacial in- 
tegral is the total momentum of the system P, which 
is the conserved quantity. By a suitable choice of the 
reference frame it can always be maid to vanish. In the 
following, we will assume that for the case of unperturbed 
system that condition is satisfied. 

Energy density p H (r;t) can be introduced in a simi- 
lar way i.e., that its spatial integral is the total system 
energy, or the Hamiltonian H, and is therefore the con- 
served quantity. It is convenient to distinguish explicitly 
contributions from kinetic and interaction parts, so that 
the microscopic expression for p H (r;t) reads [64] 



3p?,i(*) 



(2.3a) 



Pn(r;t) = p H ,kin(r;i) + p H ,int(r;i), 

p H , ki „M) = E%^,M)> 



(2.7a) 
(2.7b) 
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p H ,int(r;i) 



) Q7 -.( 



' (r,r';£) ; 
(2.7c) 



where fuv,ij( r > r '>t) ^ s a two-particle distribution func- 
tion. Its general expression depends on many factors of 
the system and can be obtained, in principle, from the 
kinetic theory. In our consideration, however, we will use 
the approximation 



/0H,int(r;t) 



CD,® 



a 7 



x * (r -■£,(*)) J^-r^i)) 



(2. 



where g"2ij(\ r ~ r 'l) i s the radial distribution function 
that depends only on the mutual distance between sites. 

Equations of motion for dynamical variables p(r;i), 
|(r; t) and /0 H (r; t) can be obtained with the use of Hamil- 
ton equations (2.3). In the reciprocal space they take the 
form of algebraic equations, since k-dependent quantities 
ip(k;t), jj(k;t) and (0 H (k;t) are spatial 3d-Fourier trans- 
forms of their counterparts in the direct space (we use 
the same notations in both direct and reciprocal spaces 
and hope there will be no confusion since we distinguish 
functions by their arguments). In particular, equations 
of motion for quantities in concern read 

p M (k;i) = ik-p(k;£), (2.9a) 
p(k;£) = ik : n(k;t), (2.9b) 
PhCM) = ik-j n (k;t)+ih(k;t), (2.9c) 

where II(k;i) is the Fourier transform of the stress ten- 
sor, and j H (k;£) and h(k;t) are related with the Fourier 
transform of the energy current (its kinetic and inter- 
action parts, respectively). For the moment, a detailed 
knowledge of their structure is not needed. 



III. THE PERTURBATION HAMILTONIAN 

In this section we are concerned about derivation of 
the perturbation Hamiltonian. Hence, it is worth to be- 
gin from considering the generation of ultrasonic waves 
and their propagation in a media. A typical type of the 
ultrasound source in action is displayed in Fig. 1. The 
sound field in front of a source is represented by means of 
lines of a constant phase, where the phase of a wave refers 
to the position of maximum or minimum in a wave. Ex- 
perimental observations tell [1, 2] that ultrasonic waves 
lose their intensity with penetration depth rapidly. It 
means that sonic amplitude in area called the far field 
is negligible. By considering a situation when the piston 
source oscillations produce the ultrasonic wave of rela- 
tively small amplitude, one may expect that the theory 
of linear response [65, 66] can be applied to describe the 
introduced perturbation almost everywhere in front of a 
piston starting from the near field and spreading forward 
to the area beyond the focus. 



Piezoceramic 
transducer 




Near field 



Region of 
confusion 



Focus 




Far field 



FIG. 1: Piston source operating with the diameter (apertu- 
re) very much greater than the acoustic wavelength in the 
medium in contact with the front of the transducer [4]. The 
sound field in front of a piston source is drawn as lines of a 
constant phase, i.e., the positions of maximum or a minimum 
in a wave. 



When the external pressure is applied to the system, 
the perturbed part H' of the total Hamiltonian H tot can 
be related to the work on changing of the deformation 
tensor ups (here indices (3, S are used to denote Cartesian 
coordinates x, y, z). In the linear approximation this 
tensor is defined as [3] 



U/3S 



1 fdl 



£ + 9ls \ 

2 \dr s drp) ' 



(3.1) 



where 1 is the time-dependent displacement due to the 
periodically applied external pressure (which for definitc- 
ness we will think is applied along the z-axis of the sys- 
tem), and therefore should be presented as 



1 — * I(t) = acxp {— i (cut + ip)} 



(3.2) 



with a — |a|, u> and ip being the amplitude, angular frequ- 
ency, and initial phase, respectively. Assuming that the 
displacement \(t) is sufficiently small (in the macroscopic 
scale) , elementary work on changing the deformation ten- 
sor can be calculated as the internal stress force times the 
displacement. Total work is then obtained by integration 
over the volume of the system: 



dV" W(r;t) = / dV^ 



d[U(r;t)} 



dM 6 



;i(i)L, (3.3) 



where W(r; t) denotes the work of internal stress forces 
in unit volume. The density of this work is the density 
of the energy of elastic deformation and can be regarded 
as the density of the perturbation Hamiltonian H'(r;t). 
For the following consideration it is convenient to rewrite 
it as 



H'(r;t) = -M(r;t)E(t), 



(3.4a) 
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where 



M(r;t) 
E{t) 



/3,<5 L ,s 

exp {— i (u)t + (p)} . 



(3.4b) 
(3.4c) 



The reason for that will be seen later. 



IV. LINEAR RESPONSE 
TO THE ACOUSTIC PERTURBATION 

In the theory of linear response [65, 66] the total 
system Hamiltonian H to t is considered to be a sum of 
its regular part, say H, and perturbation, say H' . If 
the density of the perturbed part can be written as 
H'(r;t) = —M(r;t)E(t), then in terms of the spatial 
Fourier components the time-dependent average linear 
response of a dynamical quantity A(t) produced by the 
applied field E(t) is given by 

(A(k;t)) =-j±= I ds(A(k;t-s)M(-k;0))E(s), 

(4.1) 

where fee is the Boltzmann constant, T is thermodynamic 
temperature, and the subscript is used to indicate the 
average of a perturbed quantity. An example: E is an 
electric field and M is the electric dipolc moment. Both 
M and E terms for our case are defined in the previous 
section. 

By imposing fluctuation notations, p(k; t) = p(k; t) — p, 
and taking A(k; t) = p(k; t), one obtains 



(p(k;i)) c 



-i(ajt+ip) 

k B T 



x / , °°dre^^iA ;5 a /3 (p(k;T)[n*(k;0)] /35 ), (4.2) 

P,s 

where * means complex conjugation. Integrand in the 
equation above can be transformed with the use of prop- 
erties of time-correlation functions [67-70] and equations 
of motion (2.9a) and (2.9b). The result can be written 
compactly as 



<p(k;*)>« 



m k -i(ut+tp) 

k B T 



D p (k;uj) (4.3) 



with Ofe = a/k and 

fiiy,a-y 

x {fc 2 [J^AOL +«?S%(k) +iu 3 F?J(k;u;)} , (4.4) 

where S^J(fc) is the site-site static structure factor, 

F"J(k; lo) is the upper- half Fourier transform of the site- 
site intermediate scattering function F£J(k;t) with re- 
spect to time t, and J"J(fc) is tensor of site-site static 



current correlation functions, all of them introduced in 
the standard way: 

S%{k) = (N^y^ipl'*^; Ojp^kjO)), (4.5a) 

i,j 

F«J(k:r) = (JV M JV v )-'g(^(k;0) / ^ ii (k;r)),(4.5b) 

i,j 

K~m] ps = (^)"'E<[j^( k )]«[j?.*( k )]^ 4 - 5c ) 

i,3 

and 

/•OO 

F%?(k;u) = drc— F«J(k;r). (4.6) 
Jo 

We also assume for simplicity that the wave-vector k is 
taken to be parallel to the z-axis, i.e., k = (0, 0, k). 

The applied field to the system (in our case the ap- 
plied pressure) must in general be real, so that the full 
monochromatic "force" E(i), Eq. (3.4c), should be the 
superposition (as well as the displacement (3.2)) 



- (e-^ 4 + e tuJt ) = cos (cut + v ) 



(4.7) 



and the total response function is the superposition of 
responses from each component [65, 66, 68, 69]. It is 
quite straightforward then to show that 

(p(k; i)) 0l totai 

= { e ~ iM+V) ^( fc ; + e i( " t+ ^(fc; -«)} . 



(4.8) 



In the following, we will not use the subscript "total" 
next to the symbol by assuming that the quantity in 
concern indicates already the superposition of responses 
as has been explained above. 



V. ENERGY DISSIPATION 

The response of the real system to the external mo- 
nochromatic perturbation is accompanied by the absorp- 
tion and propagation of energy The reason for that is 
because under the influence of the external perturbation 
the system changes its state [68] . The difference between 
the energy absorbed and emitted is the energy dissipa- 
tion. The energy Q(\c;u)) dissipated in the unit of time 
is the averaged over one period of monochromatic field 
the time rate of change of the system energy, and for our 
case can be written as 



Q(k;w) 



2tt 



dt (M(k;i)) ( 



dE{t) 
dt 



(5.1) 



with M and E having their meaning given by definitions 
(3.4). In terms of the linear response it produces 
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Q(k;w) 



4fc B T 



dr {c WT + e- luJT ) (M*(k; 0)M(k; r)}, (5.2) 



where exponential factors are left being not reduced as 
in (4.7) intensionally in order to show the contribution 
from each term. The expression above is calculated as 

Q(k; W ) = ^ E "»>2(W)' 

x{^7(k, W )+F / S7(k,-o;)} ) (5.3) 

which means that the energy dissipation in the linear 
response theory for a system like our is determined solely 
by the site-site intermediate scattering functions of the 
system alone. 

When the sound wave is propagated through the liquid, 
its intensity decreases with the distance. In particular, 
the existence of viscosity and thermal conductivity re- 
sults in the dissipation of energy in sound waves, and the 
sound is consequently absorbed. The decrease will occur 
according to a law e~ 2ar , and the amplitude will decrease 
as e _Qr , where the absorption coefficient a is defined in 
terms of rate of energy dissipation Q and the density of 
the energy of the sound wave Q as [36] 



a=\Q\/2cQ, 



(5.4) 



where c is the (adiabatic) velocity of sound in the liquid, 
and Q is given by Q = -^PmVqV with w being velocity 
of a fluid in the sound wave. Following the definition, 
let us introduce generalized (wave-vector and frequency 
dependent) absorption coefficient a(k;cj) by the relation 



a(k; oj) 



|Q(k;a;)| 
2w s (k)QH' 



(5.5) 



where u s (k) denotes k-dependent (adiabatic) sound ve- 
locity, and we use Q — > Q(oj) = ^p m Vq(u>)V, in which 
vq — ► vo(w) = max [— aw sin (tot + ip)] = sllo is the ve- 
locity of fluid due to the external perturbation and not 
due to the thermal motion. When Q(k;u;) and Q(w) are 
substituted into (5.5) the result is 



a(k;w) = 



Ak 2 k B Tp M Vs(k) 



E 



(5.6) 



This is the first ever expression for ultrasonic absorption 
coefficient in liquid obtained from the microscopic, i.e., 
statistical-mechanical theory. Its frequency dependence 
is determined by the product of ui 4 and site-site dynamic 
structure factors F"J(\c,±w). Also the result itself is 
independent of oscillation amplitude a as it should be. 

Since the amplitude characteristics are related by lin- 
ear relations, the exponential law of attenuation with 



generalized absorption coefficient (which also can be 
called as the generalized amplitude attenuation or gener- 
alized spatial attenuation coefficient) a(k;u>) is valid for 
any related acoustic parameter. For example, the gen- 
eralized damping time coefficient that characterizes the 
damping of the wave in time is introduced as r(k; u>) = 
[u s (k)a(k; , so that the generalized temporal at- 
tenuation coefficient becomes 0(k;w) = [r(k;w)] _ = 
i> s (k)a(k; u>), and the generalized logarithmic damping 
decrement becomes i? = &6(k;w), where & is the pe- 
riod of ultrasonic wave. 

We see, then, that the detailed analysis of expression 
for either (p(k;t)} , Q(k;cj), or a(k;w) is possible after 
evaluation of D p (k; ui) . This is considered in the following 
sections. 



VI. TCFs AND THE SITE-SITE MEMORY 
EQUATION 

Frequency/time dependence of D p (fc;w), Eq. (4.4), is 
determined by the term m'^TJF"J(k;Lo), i.e., by fluc- 
tuations in number-, charge-, and mass-densities of the 
system. A suitable way to calculate such time-correlation 
functions is based on the approach by Mori [71, 72]. Al- 
though we follow it in this section, we shall not repeat 
the very details of this method: they can be found either 
in the original source or in numerous publications after. 
But definition and handling of a set of the so-called slow 
variables of the system are given, which is required at 
least to establish notations. 

Let us consider a set of dynamical variables 

a LiH = { {p£(k;t)} , {#, L (k;i)} ,/3 H (k;t)}, (6.1) 

i.e., the set consisting of three subsets of fluctuations of 
site number-densities, longitudinal site current-densities, 
and fluctuation of energy density, respectively. In the fol- 
lowing, it will be desirable to distinguish between these 
three subsets by assigning them consecutive labels (1), 
(2), and (3). From the time-/space-inversion symme- 
try properties of elements of the set 2l L>H one can see 
that J? L (k) is orthogonal to both /3"(k) and /3 H (k), while 
/3"(k) and p H (k) are not mutually orthogonal. For the 
sake of simplicity of future considerations it is advan- 
tageous to replace p H (k) with its renormalized variable 
p E (k) defined as 



p E (k) = p H (k) - e 



x / 5«(k)[S- 1 (fc)]J(^*(k) / 5 H (k)), (6.2) 

so that 

(pr(k)/5 E (k))=0. (6.3) 
We will denote the new set of dynamical variables as 

2l L;E = { {p«(k;t)} ,{j° L (k;t)} ,p E (k;i)}- (6.4) 
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The microscopic equations of motion for the elements of 
the set 2l LjE read 

3°(k;t) = ifc#, L (k;t), (6.5a) 
i° L (k;t) = z/cIT« zz (k;t), (6.5b) 
/3 E (k;i) = ifcj B>L (k;t) + ft(k;t). (6.5c) 

where H^ zz (k;t) and j EjL (k;i) are defined by 

II M (k;t) = ]T™«n^(k;t), (6.6a) 

jB,L(k) = j H ,L(k)- X! (*W*#,lO<) 

/j.,i>;a,7 

x[S- 1 (fc)]^<p2; , *(k)^(k)>. (6.6b) 



One can see now that arbitrary elements from different 
subsets of 2t L . E are initially orthogonal to each other, 
while elements within a same subset are not. It implies 
that the matrix of initial TCFs has block-diagonal struc- 
ture 

/C (1)(1) (fe) \ 

C(fc) = C (2)(2) (fc) , (6.7) 

\° C (3)( 3)(fc)/ 

where are zero-matrices, and C^^(k) are square sub- 
matrices constructed exclusively with the use of elements 
from (i)th subset: 



a 7 
(1)(1) 

(X V 

a. 7 



[C(fc)](i)# = <^'*(k)p°(k)>, (6.8a) 



[C(fc)]( 2)( i) = 07;*(k)i« L (k)), (6.8b) 



[C(fc)](3)(3) = (p E (k)p E (k)). 



(6.8c) 



Important property of the block-diagonal matrix is that 
its inverse has similar block-diagonal structure: 



/C^ (1) (*0 



c 1 

*-(2)(2) 







(k) 



(6.9) 



^-(3K3)( fc ) y 



That feature will be used later. Meanwhile, one has to 
note that in the case t ^ matrix of TCFs is not nec- 
essary block-diagonal. For simplicity of notations, we 
will denote elements of the entire set 21 l ,e at time t by 
A\ (k; t) , where A is a composite label of the element of 
the set. Vector- row formed by these elements and vector- 
column formed by their complex conjugates will be de- 
noted, respectively, by A(k;i) and A*(k;i). Therefore 
the second-rank tensor of TCFs C(fc; t) constructed from 
elements of the set 21 l ,e, and its components read, respec- 
tively, 



C(M) 
[C(M)] AlAa 



(A*(k;0)A(k;t)>, (6.10a) 
(A* X2 (k;0)A Xl (k;t)). (6.10b) 



If one has to be more specific, it may be necessary to 
assign indices A their values which can be fi, a, and the 
subset number (i). For example: 



[C(k;t)} (2)(i) 



(p3'*(k)i«(k;t)), 



(6.11a) 



[C(fe;t)](i)(3) = <&(k)p°(k;t)>, etc. (6.11b) 

n 

Thus the problem of evaluation of a(k;u) amounts to 
finding out the solution for time development of TCFs 
C(i)(i)(M). 

By introducing the projection operator whose ac- 
tion on arbitrary variable is described as 

&...= ^(^(k)...}^- 1 ^)]^^^^^, (6.12) 

and its complementary projection operator =2 = 1 — & ', 
one can obtain the so-called memory equation for TCFs 
C XlX2 (k;t) = [C(fc;i)] AlAa in the form [71-74] 

C XlX2 (k;t)=^iQ XlX ,(k)C x > X2 (k;t) 

A' 

-J2 [ dTK XlX ,(k;T)Cy X2 (k;t-r), (6.13) 

A' J ° 

where i£l XlX2 (k) = [iil(k)] x x are elements of the so-cal- 
led frequency matrix given as 

*0 AlA2 (k) = Y,(A X , (k)iJSM Al (k)) [C-^fe)] v , (6.14) 



K XlX . 2 (k;t) = [K(k;t)] x x are elements of the matrix of 
the first order memory kernels introduced as 

K XlX2 (k; t) = J2(fx> (k)fA, (k; t)) [C" 1 ^)] x ,, 2 , (6.15) 



and fA(k; t) are the so-called fluctuating forces defined as 
f A (k;t) = exp{^J^}fA(k;0) 



= c 



A x (k)-J2^ xx ,(k)A x ,(k) 



.(6.16) 



In equations above, _Sf is the Liouville operator of the 
system. Explicit structure of all nonzero elements of ma- 
trices (6.14)-(6.16) is listed in Appendix A. 

Since in our case all dynamical variables are related 
to specific interacting sites, we shall call the correspond- 
ing memory-equations (and subsequently all quantities 
involved) the site-site memory-equations. Study of dy- 
namic processes in molecular liquids with the use of the 
site-site memory equations has been under development 
in last decade in works by Hirata and coworkers (see, 
e.g., Refs. [55-63] and references therein). The present 
treatment is, however, substantially different, and the 
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difference appears in the choice of the set of slow vari- 
ables: fluctuations of the total energy density arc in- 
cluded into the set for the first time. As we shall se later, 
it is essential in order to reproduce the expression for the 
so-called classical coefficient of ultrasonic absorption in 
liquids. Solving equation (6.13) may give (in principle) 
the answer about the structure of intermediate scattering 
functions F"J(k;t) through the finding C(i)(i)(fc; t). Its 
solution in the frequency domain is written briefly as 

C ( i) (1 )(fc;w) = Y- 1 (fc; W )C (1)(1) (A : ;0), (6.17) 
where C(i)(i)(fc;0) = C(i)(i)(fc) is matrix of initial valu- 



es, and V~ 1 (k;u) can be regarded as a sort of propaga- 
tor since all frequency/time dependence of the solution 
is determined by this quantity. Required steps to find 
the formal solution to the site-site memory equation as 
well as the solution itself as a function of frequency terms 
and memory kernels is displayed in Appendix B. Gen- 
eral expression for Y(fc;o>) is rather lenghty, but if time 
is long and wave-vector is sufficiently small (in the case 
of propagation of ultrasonic wave both these conditions 
are satisfied), cross memory terms can be neglected lead- 
ing to 



Y(fc;w) 



K(i)(i)(fc;w) - \iw 

-ifi(i)( 2 )(fc) |k ( 2)(2)(A;;w) - \iu - ifl {2){3) (k) K ( 3 )(3) (fc; w) - \iw «f2 (3)(2) (fc) [ iil ( ., uil (k). (6.18) 

I 



One can also replace in remaining expressions for mem- 
ory kernels the evolution operator with complementary 
projection by an ordinary one, and at the same time eval- 
uate the results at the leading order k 2 . Such operation, 
however, cannot be justified for the general case, and it 
makes sense to talk about the final solution to the me- 
mory equation if one specifies the model for the memory 
kernels. Since ultrasonic frequency is sufficiently low, the 
dynamics of the system should be considered in the long- 
time limit, i.e., in which the time scale is large enough for 
fast relaxation processes to be completed. Memory kernel 
of the ordinary memory equation in that case is usually 
constructed by the so-called mode-coupling approxima- 
tion [73, 74]. In works by Chong et al. [75, 76] the 
conventional mode-coupling theory has been extended to 
the case of molecular liquids based on the interaction-site 
model. One shall follow this procedure for the numerical 
evaluation. Albeit, some analytical analysis is possible 
for the case of the one-component simple liquid in the 
hydrodynamic limit. 



VII. INITIAL VALUES AND FREQUENCY 
TERMS 

In order to be solved, either analytically or numeri- 
cally, site-site memory equation (6.13) requires informa- 
tion about initial value. Since we are interested in site- 
site intermediate scattering functions, the initial values 
of our concern are site-site static structure factors. They 
can be determined within the frame of the integral equa- 
tion theory of liquids called RISM, or one of its general- 
izations (e.g., extended RISM [48-51], etc.). It predicts 
static structure of molecular fluids via the calculation 
of site-site pair correlation functions. This method has 



been extensively used and proved to be the powerful tool 
in the microscopic description of equilibrium qualities of 
the system. The main equation for mixture can be writ- 
ten in the reciprocal space in matrix notations as 

php = wcw + ujcphp, (7.1) 

where p is diagonal matrix consisting of number densities 
of each molecular species, while h = h(fc), c = c(fc) and 
lj = u>(k) are matrices of the Fourier transform of site- 
site total, direct and intra-molecular correlation func- 
tions, respectively. Equations (7.1) are solved with the 
closure specified. Typical closures are hyper-netted chain 
(HNC), mean-spherical approximation (MSA), Percus- 
Yevick (PY), etc. [51]. Then the matrix of site-site struc- 
ture factors is given by 

S(k) = p- 1 uj{k) + ph(k). (7.2) 

The static site-site current correlation function can be 
treated analytically. The expression for arbitrary shape 
of the molecule has been given recently [60, 61] and reads 

[j L (fc)n - ^jo(fc^) 

x[Sr a ] T ■ [(Trl- 1 )!-^ 1 ] • tfr* 
- ^yh{k^)[5v^5v^Y ■ I 1 ■ [Sr^nSr-y], (7.3) 

where M is the total mass of the molecule, l ai is the 
distance between sites, 5r a is the vector pointed from 
the center of mass to the site a, I is matrix of iner- 
tia moments of the molecule, I has the same meaning of 
the diagonal unit matrix; finally, j and j 2 are spheri- 
cal Bcsscl functions of the first kind of zero and second 
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order, respectively, and symbol tt is used to denote the 
outer vector product. 

Expressions (7.2) and (7.3) are crucial in order to cal- 
culate the solution (6.17), linear response terms (4.4), 
and frequency terms (6.14). Actually, calculation of fre- 
quency terms i£la)( 2 ){k) and i£lr 2 )(i){k) is rather triv- 
ial, while calculation of i£l( 2 )(3)(k) and if2( 3 )( 2 )(fc) may 
be difficult. The difficulty is twofold: one is caused by 
b^(k), and the other by (/5*(k)/5 E (k)). In the former case 
one has to deal with three-body correlation functions, 
i.e., 03*(k)II« zz (k)) and (/}*(k)^(k)), while in the lat- 
ter case it is essentially the problem of calculation of 
four-body correlation function. It follows from the defini- 
tion of the microscopic energy density of the system, eqs. 
(2.7) and (2.8), where interaction part involves the two- 
body distribution function. There is no reliable mecha- 
nism in statistical mechanics which is capable to handle 
this accurately for the general case. Experiments also do 
not reveal many-body correlations directly (unlike pair 
correlations) [77]. If one does not wish to involve com- 
puter simulations, all that leaves quite limited space to 
maneuver and narrows numbers of choices to just one: 
approximations. Particular approximations may be case 
dependent, but fortunately there is one that may work for 
arbitrary system. It is given by hydrodynamic limit and 
for most cases can be calculated nearly exactly. Therefore 
approximation of b"(k) and (/5*(k)/5 E (k)) by their hydro- 
dynamic values sounds naturally in view of the slowness 
of ultrasonic processes and is considered as an effective 
choice for numerical evaluations. 



VIII. RELATION WITH HYDRODYNAMICS 

In order to demonstrate the robustness of our theory, 
here we examine its hydrodynamic limit. For that pur- 
pose we confine ourselves to the case of one-component 
simple liquid. The expression for the so-called classi- 
cal ultrasonic absorption coefficient in liquids a ch derived 
within the framework of a continuous media theory is well 
known [1, 2, 36]: 



2p M v% 



3^s 



Here v s is (zero frequency) adiabatic sound velocity in 
the liquid, ?7 S and n v are phcnomcnological coefficients of 
shear and bulk viscosities, respectively, A T is phcnomcno- 
logical coefficient of its thermal conductivity, and Cy and 
Cp are specific heat capacities per unit mass of the liquid 
at constant volume and pressure, respectively. 

The starting point to relate our theory and hydrody- 
namics would be the formal solution to the equation for 
intermediate scattering function specified by eqs. (6.17) 
and (6.18) in which Y(k;uj) for one-component simple 
liquid reduces to 



Y{k-uj) 



iQ,\ 2 (k) i$l 2 i{k) 



-iuj + K 22 (k;uj) 



in 23 (k) in 32 (k) 



1.2) 



-iu + K 33 (k;uj) 
Expressions for the frequency terms in that case read 



iSli 2 (k) — ik, 
iil 2 i(k) = ik 

ifl 23 (k) 

ifl 32 (k) 



k B T 
mS{k) ' 
ik b{k) 
m{p*(k)p E (k))- 
ikb*{k) 



(8.3a) 
(8.3b) 

(8.3c) 

(8.3d) 



Nk B T ' 

where S(k) is the static structure factor of the liquid, and 



b(k) = (&(k)n„(k)) 



k B T 
S(k) 



</3*(k)p(k)). (8.4) 



In the hydrodynamic, i.e., k — ► limit 

S(k)\ k ^ = k B Tpx T , 



(8.5) 



b(k)\ k ^ Q = k B T V — , (8.6) 

Xt 

(p E (k)p E (k))| fe ^ - k B T 2 Nc v , (8.7) 

where Xs,t and (3 P are adiabatic/isothermal compressibil- 
ity and isobaric thermal expansion coefficients, respec- 
tively, defined in thermodynamics as 



1_ fdv\ 



0* 



1 fdV 



V \ dT ) p ' 



(8.8) 
(8.9) 



and c PiV is specific heat capacity per particle at con- 
stant pressure/ volume. To complete the treatment of fre- 
quency terms one may also need the thermodynamic de- 
finition of adiabatic/isothermal sound velocity v SjT , i.e., 



f dP 



\9pM/gx PmXs.t 



(8.10) 



and thermodynamic expression that relates heat capaci- 
ties at constant pressure and volume, i.e., 



Cp C v 



T0- 



PXt 



m) 



With all that in mind, the corresponding products of 
frequency terms result in 



in 12 {k)ifl 21 (k)\ k ^o = -k 2 ^, (8.12a) 

7 

in 23 (k)m 32 (k)\ k ^ = -k 2 v 2 s (i - , (8.12b) 



-iuj 



where 7 = c P /c v = Xt/Xs- 
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Remaining memory kernels K 2 2{k;w) and K 2 2(k;Lj) 
should be treated as follows. It is known [73, 74] that ge- 
neralized longitudinal viscosity, shear viscosity, and ther- 
mal conductivity coefficients, respectively, can be written 
in terms of familiar TCFs as 



F(k;u,)| hyd = F(k; 0)| hyd 



r/ L (k;t) oc 



1 



(n: z (k;0)n zz (k;*)}, (8.13a) 



k B TV 

%(k;t) K fc^ n -( k ;°) n -( k ;*)>< ( 8 - 13b ) 

A T (k;t) cx ^l^(j* L (k;0)j BlL (k;i)>. (8.13c) 

These expressions usually are complemented by the defi- 
nition of the generalized bulk viscosity coefficient ?7v(k; t) 

as 



ry L (k;t) = - Vs (k; t) + 7] v (k; t). 



(8.13d) 



The proportionality sign in (8.13a)-(8.13c) can be re- 
placed by the equality sign after specifying the conven- 
tion used to identify the scalar product of two dynamical 
quantities, say A and B, with their time-correlation func- 
tion. The standard correspondence (A, B) = (A*B) — 
(A*)(B) is due to Mori [71] and has to be used if both 
A and B have nonzero averages. This is precisely hap- 
pening in the case for r] L (k;t), since [73] (U zz (k))\k^o = 
(2n) 3 5(k)pkBT. Changing from the continuous to the 
discrete reciprocal space, the correct definition for f] L (k; t) 
is obtained as 



?7 L (k;t)| fc ^o = 



1 



k B TV 

AS£t 



xe^ t [n zz (k^0)-PV}), (8.14) 



while both r/ s (k;t) and A T (k;i) do not require correction 
since (Il zx (k))\ k ^o = by definition and (j EiL (k))| fc ^ = 
by the convention used in this paper (that was discussed 
in Section II). Relations of generalized transport coeffi- 
cients 7(k;t) = {r] L (k;t),r] s (k;t),f] v (k;t),X T (k;t)} with 
ordinary hydrodynamic longitudinal, shear, and bulk vis- 
cosities, and thermal conductivity T = {i] L , ?7 S , rj Wl A T }, 
respectively, read 
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f-OO 

= dt limT(k;£), 
Jo 



3»7s + Vv 



(8.15a) 
(8.15b) 



Hence, in the hydrodynamic limit memory kernels are 
identified as 

K 2 2(k;u))\ k ^ ,^o = k 2 v h , (8.16a) 
k 33 (k;uj)\ k ^ Q ^ Q = k 2 D Tl , (8.16b) 

where v L — rj L /p M is the longitudinal kinematic viscosity, 
and D T = A T / pc P is the thermal diffusivity. 

When expressions (8.12) and (8.16) are substituted 
into (8.2), the solution for F(k;u) in the hydrodynamic 
limit is obtained as 



-, -l 



-iui + 



k 2 v 2 s h 



■ , 12 , * 2 «s (1-1/7) 

— lu + k L> T 7 . 



(8.17) 



where .F(k;0)|h y d = S(k)\k^o an d is given by Eq. (8.5). 
By using that to calculate Q(k;cj)|h y d and saving the 
only leading order contributions, the generalized ultra- 
sonic absorption coefficient in the hydrodynamic limit is 
obtained as 



a(k; w)|h y d = 



2p M v 3 s 



1 

Cv 



1 



(8.18) 



Taking into account the difference in definitions of Cp, 
C'v and c P , c v , one can see that it completely coincides 
with a C L given by Eq. (8.1), which means that our theory 
is hydrodynamically consistent. And the importance of 
fluctuations of the total energy density p H (k;t) in the 
set of slow variables 2l LjH now is revealed: it would be 
impossible to obtain the classical ultrasonic absorption 
coefficient with contribution from thermal processes if 
p H (k; t) is not considered explicitly. 



IX. SUMMARY 

In this work we presented a statistical-mechanical the- 
ory for treatment of such ultrasonic processes as propa- 
gation and absorption of ultrasound in molecular liquids. 
In particular, we demonstrated the way to calculate gen- 
eralized, i.e., wave- vector and frequency dependent ul- 
trasonic absorption coefficient a(k;u>). The suggested 
description is a combination of the linear response the- 
ory by Kubo [65, 66] and the memory equation formal- 
ism by Mori [71, 72] extended to the case of molecular 
liquids based on the interaction site model [51]. Our re- 
sult is the first one obtained from the microscopic theory, 
and reflects all the general peculiarities of ultrasonic pro- 
pagation and absorption in liquids known from the phe- 
nomenological description. According to classification by 
Dukhin and Goetz [6], these peculiarities can be sorted 
in several different mechanisms. 

The viscous mechanism is hydrodynamic in nature. It 
is related to the shear waves generated by particles or 
groups of particles oscillating in the acoustic pressure 
field. The shear waves appear due to the difference in 
densities in the vicinity of particles and next nearest 
medium. For example: in terms of Fig. 1, lines of con- 
stant phase may be thought as positions of maximum 
compression, while places in between them - as posi- 
tions of dilations. As a result, the liquid layers in the 
particle vicinity slide relative to each other, and the non- 
stationary sliding motion of the liquid near the particle 
is referred to as the shear wave. The viscous mecha- 
nism is considered to be the most important for acoustics 
since it causes losses of acoustic energy due to the shear 



11 



friction. In our work, viscosity related contribution into 
absorption coefficients is controlled mostly by the mem- 
ory kernel (A4e), whose hydrodynamic limit expression is 
proportional to the standard time-autocorrelation func- 
tion of the stress tensor (8.14) and is therefore associated 
with generalized longitudinal viscosity. 

The thermal mechanism is thermodynamic in nature 
and is related to the temperature gradients in the vici- 
nity of particles and next nearest medium. The temper- 
ature gradients are due to the thermodynamic coupling 
between pressure and temperature. There is variety of 
experiments telling that dissipation of acoustic energy 
caused by thermal losses is much smaller than one caused 
by the viscous mechanism if it is the case of liquids and 
solutions, but not liquid metals [1, 2]. However it may 
be dominant for colloidal systems with soft particles in- 
cluding emulsion and/or latex droplets [6]. Although the 
latter systems are not a subject of our consideration, the 
thermal mechanism of ultrasonic losses in our description 
is included explicitly. One of consequences of this inclu- 
sion is the presence of thermal conductivity coefficient of 
a liquid in the hydrodynamic limiting expression of its 
ultrasonic absorption coefficient. The required condition 
for this is the consideration of fluctuations of the total 
energy density (2.7) as one of parameters of abbreviated 
description forming the set (6.1). It is worth to note that 
in the limit of small wave- numbers two major compo- 
nents of ultrasonic absorption - viscous and thermal - 
are additive, eqs. (8.1) or (8.18), while in general they 
cannot be separated from each other. The concept of a 
combined mechanistic in nature viscous and thermody- 
namic in nature thermal treatment of ultrasonic absorp- 
tion process in a liquid can be realized only in microscopic 
description of the problem as it was demonstrated in our 
paper. 

The electrokinetic mechanism describes interaction of 
ultrasound with the double layer of particles/molecu- 
les. Oscillation of the charged particle in the acoustic 
field leads to generation of an alternating electric field, 
and consequently to an alternating electric current. This 
mechanism is used in electro-acoustic measurements. In 



the interaction site model of molecular liquids individual 
sites are usually atoms or atomic groups of the molecule 
and therefore often have an electric charge. Hence, in our 
theory most of the types of electrostatic interaction are 
included through the system Hamiltonian, and most of 
the types of interaction of ultrasound with double layer 
are taken into account through the formalism of the linear 
response. In particular, the response of charge-density 
fluctuations is described by eq. (4.8) and provides the 
basis for calculation of ultrasonic vibration potential gen- 
erated in the system. 

The short wave-number limit considered at the end 
of our paper demonstrated hydrodynamic consistency 
of our calculations in a sense that we were able to ac- 
curately reproduce the result for ultrasonic absorption 
coefficient known from continuous media theories and 
called therefore the classical absorption coefficient. On 
the other hand, use of hydrodynamic limit for some fre- 
quency terms and the memory kernels, as described in 
Sections VII and VIII, is rather a necessity steamed from 
the fact that the requirement of analytical handling of 
many-body correlation functions at finite wave-numbers 
represents a considerable challenge. For all other quanti- 
ties one shall be able to use their microscopic expressions. 

We expect our theory of ultrasonic absorption to work 
for most of the types of molecular liquids, including sim- 
ple, non-polar, polar and ionic. The validity of this ex- 
pectation has to be tested by applying the presented for- 
malism to some real systems and evaluating their ultra- 
sonic absorption coefficients numerically. It could very 
well be a subject for another paper. 
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APPENDIX A: MATRIX ELEMENTS OF THE SITE-SITE MEMORY EQUATION 

Nonzero elements of the site-site frequency matrix ifl(k) read 

v' ,7' v' \-y f 



[iSl(k)] |) ( T) = £ <#,''*(k)i« L (k)} [C- X (fc)] (t) ( T) ]T (N„,N v )-i(ji: : :(k)j« h Q<.)) [S-\k)] J 2 , (Alb) 

[i«(fc)] (2)(3) - (p*(k)i« L (k)) [C-^fc)] (3)(3) =lk 6«(fc)/(p*(k) / 3 E (k)) ) (Ale) 

[i«(fc)] ( 3 )(Z) = E M/W) [C-'W] (fe =ik E {N vl N v )-h^*{k) [K\k)]t>v > ( Ald ) 
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where 



6»(fc) - (p*(k)n«, z (k)> - E ( ^ ) " (pS(k)^(k)> [s- 1 ^)]™ [j L (k)]° 



at, 



172 

I L_L,, '" / JM 1/ 2 



(A2) 



"ii"2;7ii72 

Nonzero elements of the site-site vector of fluctuating forces f(k) at initial time t = are then obtained in the form 



[f(k)] (!) = l k 

ft 



[f(k)] (2) = 
A* 



[f(k)] (3) = ik 



n^(k)- E f^)V(*)]# ^(^^^(k) 



bg(fe)fe(k) 

(^(k)p E (k)) 



jE , L (k)- 2 {N v ,N v )-hi'\k)[}-\k)Y2j2 t ^) 



u,u' ;7,7' 



+ /i(k) 



=iAj° L (k), (A3a) 

^n^(k), (A3b) 

EifcJ EiL (k). (A3c) 



In terms of elements of a set of dynamical variables 2l LjE , Eq. (6.4), site-site frequency matrix itt(k), eqs. (Al), and 
site-site vector of fluctuating forces f(k), eqs. (A3), elements of the matrix of the first order site-site memory kernels 
K(k;t) are expresses as 



K(M)] 


a 7 

(i)(i) 

j-L V 


= fc 2 E WW* <^;;(k)e^« t (k)> [S-^Jfe)]^ , 

I/' ,7' 


(A4a) 


K(k;t)} 


a 7 
(1)(2) 


= fc 2 ]T <nj; ; ^(k)e«^*i- t (k)> [J^W]^, 

1/' ,7' 


(A4b) 


K(k;t)} 


(1)(3) 


= fc 2 0- L (k)e^ t J« L (k))/( / 5 E (k)p E (k)), 


(A4c) 


K(k;t)} 


a 7 
(2)(1) 

fl V 


= fc 2 E (AW" 1 <#;* (k)e^n^(k)) [s- 1 ^)]^ , 

I/' ,7' 


(A4d) 


K(k;t)} 


a 7 
(2)(2) 

j-L V 


= fc 2 E {N v N v ,y h <n^(k)e^'n«^k)) [j- 1 ^)]^ , 

v' ,7' 


(A4e) 


K(k;t)} 


(2) (3) 
fJ- 


= fc 2 <j E * L (k)e^*n« zz (k)) / (p* (k)p E (k)> , 


(A4f) 


K(k;t)} 


(3)(1) 

V 


-fc 2 E (^)-* <^;;(k)e^ lt (k)> [S- 1 ^)]^ , 

1/' ,7' 


(A4g) 


K(k;t)} 


(3) (2) 

V 


= fc 2 E (aw- (n3:;L(k)e^\ L (k)) [j^w]:;: , 

v' ,7' 


(A4h) 


K(k;t)} 


(3)(3) 


= fc 2 (j E * L (k)e^ t J E , L (k))/( / } E (k)p E (k)). 


(A4i) 



APPENDIX B: FORMAL SOLUTION 
TO THE MEMORY EQUATION 

We are interested in finding an expression for intermedi- 
ate scattering functions F"J(k;t), i.e., C(i)(!)(fc;i) if it is 
in matrix notations. It is convenient to do calculations 
in the reciprocal space, after applying to the memory 
equation the upper-half Fourier transform. There will 
we three type of equations for C^^(k;u)): 

-iwC(i )( i)(A;;w) = C(i )( i)(fc;0) 



+zfl(i)(2)(fc)C( 2 )(i)(fc;w) 
-K(i)(i)(fe;w)C(i) ( i)(A;;a;) 
-K(i)(2)(fc;w)C(2)(i)(A;;w) 
-K(i)(3)(fe;w)C( 3)( i ) (A;;a;), 



(Bla) 



-iwC( 2 )(i)(fc; w) = ifi( 2 )(i)(A;)C(i)(i)(fc;w) 
+^(2)(3)(fc)C( 3 )(i)(fc;w) 
-K( 2 )(i)(fc;w)C(i)(i)(A;;a;) 
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-K(2)(2)(fc;w)C( 2 )(i)(fc;u>) 
-K(2)(3)0;w)C( 3 )(i)(fc;u>), 



By substituning into (Bla) expression for C(2)(i)(fc;w) 
(Bib) from (Bib) and expression for C^^(k;u)) from (Blc), 
the formal solution for C(i)(i)(fc;w) can be written as 



C(i )( i)(fc;a;)= Y(fc;co) C (1)(1) (fc; 0), (B2) 



-iwC( 3 )(i)(fc;w) = in(3)(2)(fc)C( 2 )(i)(fc;w) 
-K( 3 )(i)(fc;w)C (1 )(i)(fc;w) 
-K ( 3 )(2 )(fc;u>)C ( 2)(i)(fc;w) 

— K( 3 )(3)(fc; w)C( 3 )(i)(fc; w). (Blc) where matrix Y(fc; lo) has the form 
I 



-l _ 



Y(fc;w) = K (1)(1) (fc;w) - liw - K (1)(3) (fc; w) K (3)(3) (fc;u) - liw K (3)(1) (fc; u) 

iO(i) (2 )(fc) - K (1 )( 2 )(fc;w) - K( 1)( 3)(fe;w) K (3)(3) (fc;w) - \iu> ifl {3){2) (k) - K (3)(2 ) (fc; w) 



x <^ K(2)(2) (k; uj) - liw - ifi( 2 )( 3 ) (&) - K (2 )( 3 ) (fc; w) K( 3 )(3) {k; w) - liw in (3)(2 ) (fc) - K( 3 )(2) w) 



x < *n (2 )(i)(fc) - K( 2 )(i)(fc;w) - ifi( 2 )(3)(fc) - K( 2 )(3)(fc;w) K(3)( 3 )(fc;a;) - liw K (3)(1) (fc; w) 



(B3) 
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